All Six Planets Known to Orbit Kepler- 11 Have Low Densities 
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ABSTRACT 

The Kepler- 1 1 planetary system contains six transiting planets ranging in size from 
1.8 to 4.2 times the radius of Earth. Five of these planets orbit in a tightly-packed con- 
figuration with periods between 10 and 47 days. We perform a dynamical analysis 
of the system based upon transit timing variations observed in more than three years 
of Kepler photometric data. Stellar parameters are derived using a combination of 
spectral classification and constraints on the star's density derived from transit profiles 
together with planetary eccentricity vectors provided by our dynamical study. Com- 
bining masses of the planets relative to the star from our dynamical study and radii of 
the planets relative to the star from transit depths together with deduced stellar proper- 
ties yields measurements of the radii of all six planets, masses of the five inner planets, 
and an upper bound to the mass of the outermost planet, whose orbital period is 118 
days. We find mass-radius combinations for all six planets that imply that substantial 
fractions of their volumes are occupied by constituents that are less dense than rock. 
The Kepler- 11 system contains the lowest mass exoplanets for which both mass and 
radius have been measured. 

1. Introduction 

Within our Solar System, Earth and smaller bodies are primarily rocky (or, far from the Sun, 
mixtures of rock and ices), whereas the cosmically-abundant low-density constituents H 2 and He 
dominate the volume in Uranus/Neptune and larger bodies. There are no local examples of bodies 
intermediate in size or mass between Earth (1 R ffi , 1 M©) and Uranus/Neptune, both of which are 
larger than 3.8 R© and more massive than 14 M©. However, observations of extrasolar planets are 
now filling this gap in our knowledge of the mass-radius relationship of planetary bodies. 
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To date, the only accurate radius measurements for exoplanets have been provided by planets 
observed to transit across the disk of their star. The fractional depth of the transit provides 
a direct measure for the ratio of the radius of the planet to that of its star. The star's radius 
is estimated using spectroscopic classification, in some cases augmented by other techniques. 
Doppler measurements of the variation of a star's radial velocity have been used to compute mass 
estimates for almost two hundred transiting giant planets, as well as for the f irst three sub-Ur anus 



exop lanets for whi ch both radii and m asses were determ ined. GJ 1214 b (ICharbonneau et al 



2009), CoRoT-7 b (lOueloz et al 



20091), and Kepler-10 b (IBatalha et all 



201 Analysis of transit 



timing variations (TTVs) resulting from mutual planetary perturbations provide d dynamical 



estima tes of the masses of the five innermost known planets orbiting Kepler- 11 (ILissauer et al 



201 lal) . more than doubling the number of exoplanets less massive than Uranus with both size 
and mass measurements. Precise mass estimates have subsequently been obtained for sever al 



more sub-Uranu s mass planets , in three cases by using RV, 55 Can cre e (IWinn et al 



Endletal 



20121) . Kepler-20 b (Gautier et al 



using TTVs, Kepl er-36 b.c (ICarter et al. 



one, Kepler- 18b (iCochran et al 



20121) . and GJ 3470 b (Bonfils et al. 



20121) . and Kepler-30 b (lSanchis-Ojeda et al 



20111: 



20 



2Kthree 



2012), and 



20111) . using a combination of RV and TTV data. Less precise 



estimates for the masses of dozens of Kepler plan ets and planet candidate s, many of which are in 



this mass range, have been derived from TTVs by IWu & Lithwickl (|2012|) . 



Lissauer et al. 



(|2011a|) estimated the masses of the five planets Kepler- 1 1 b-f using only the 



first 16 months of Kepler data. Similar mass constraints on these plane ts, as well as an upper lim it 



of 30 M m on the mass of t he outer planet Kepler- 1 1 g, were obtained by 



Migaszewski et al. 



Migaszewski et al. 



(2012). 



(|2012l) analyzed the same Ql-Q6jdata using a photodynamical model, which 



'The Kepler spacecraft rotates four times per orbit to keep the sunshade and solar panels ori- 
ented properly. Targets are imaged on different parts of the focal plane during different orienta- 
tions. The Kepler orbital period is ~372 days, and the data are grouped according to the "quarter" 
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adjusted planetary parameters (size, orbital elements, masses) to minimize the residuals of a fit of 
a model lightcurve that accounts for mutual planetary interactions to the measured lightcurve. 

We report herein more precise estimates of the masses of the six Kepler- 1 1 planets derived 
from TTV measurements that incorporate 40 months of Kepler photometric time series data. 
In Section 2, we present our estimates of transit times; detailed descriptions of the three 
independent techniques used to compute these times are given in Appendix A. Our dynamical 
analysis of the Kepler- 1 1 system based upon these transit times is presented in Section 3, with 
additional information provided in Appendix B. In Section 4, we combine estimates of stellar 
density obtained using transit profiles and the dynamical measurement of planetary eccentricities 
presented in Section 3 together with analyses of high-resolution spectra taken at the Keck 
I telescope to provide refined parameters for the star Kepler- 1 1 . We tabulate the properties 
of Kepler- 11 's six known planets that are derived by combining lightcurve analysis with our 
dynamical results and stellar parameters in Section 5, wherein we also discuss implications of 
these results for planetary compositions. We conclude the paper with a summary of our principal 
results. 

2. Measurement of Transit Times from Kepler Photometric Time Series 

Variations in the brightness of Kepler- 11 have been monitored with an effective duty cycle 
exceding 90% starting at barycentric Julian date (BJD) 2454964.512, with all data returned to 
Earth at a cadence of 29.426 minutes (long cadence, LC); data have also been returned at a 
cadence of 58.85 seconds (short cadence, SC) since BJD 2455093.216. Our analysis uses short 
cadence data where available, augmented by the long cadence dataset primarily during the epoch 

year during which observations were made. The data on Kepler- 11 taken prior to Kepler's first 
"roll" are referred to as Ql. Subsequent quarters are numbered sequentially: Q2, Q3, ... 
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prior to BJD 2455093.216, for which no SC data were returned to Earth. We obtained these data 
from the publicly-accessible MAST archive at http://arch ive.stsci. edu /kepler/ . 

As measurement of transit times (TTs) requires a complicated analysis of often noisy data, 
authors Jason Rowe (J.R.), Eric Agol (EA.) and Donald Short (D.S.) performed independent 
measurements of TTs using techniques described in Appendix A. Figured] shows the deviations of 
all three sets of observed transit times, O, relative to time from a linear ephemeris fit, Ci, through 
Q14 Kepler data. Here and throughout we base our timeline for transit data from JD-2,454,900. 

As evident in Figured] each set of TT measurements contain several outliers. These outliers 
are unlikely to be correct, and may be due to overlapping transits, star spots, or uncertain fits to 
the lightcurve. Trying to fit these outlier TTs would degrade our dynamical studies. Therefore, we 
remove points where only one of the methods yields a TT whose uncertainty is less than 2.5 times 
as large as the median TT uncertainty computed by that method for the planet in question. We 
then use the three sets of measured TTs to filter out unreliable measurements as follows: If two or 
three sets of measurements are available for a specific transit and each of the la uncertainty ranges 
overlap with at least one of the other ranges, then each of the points are used. If there is only a 
single measurement, or if there is no overlap of la uncertainty ranges, then all measurements 
of this transit are discarded. If three measurements are available, and two overlap but the third 
does not overlap with either, then the data are examined by eye to determine whether TTs of the 
overlapping pair are sufficiently robust to include. This culling procedure removed fewer than 
8% of detected transit times from each dataset, with the most points discarded from Kepler- 11 b, 
whose transits are the most numerous and have the lowest signal-to-noise ratio (S/N). For planet 
b, we removed 17 of the 103 TTs measured by E.A., 9 of the 111 TTs measured by J.R., and 13 of 
the 90 TTs measured by D.S. Our approach is conservative in the sense that the data set used for 
our dynamical studies presented in Section 3 consists only of transit times that are corroborated 
by at least one alternative method. 
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Fig. 1. — Transit timing variations for Kepler- 11 's six known planets, using short cadence data 
when available, supplemented by long cadence data prior to £(JD-2,454,900) = 193 days, where 
short cadence data were not sent to Earth. The TTs measured by E.A. are displayed as green 
open triangles, those from J.R. as blue open circles, and those calculated by D.S. as red open 
squares, with their respective methods described in Appendix A. The sets of data points are largely 
consistent. The observed transit times, O, are displayed as deviates from times, Cu that were 
calculated using a linear fit to each set of transit data, i.e., a fit that assumes strictly periodic orbits. 
All measured TTs are displayed, apart from one outlier for Kepler- 1 1 d that deviated from both of 
the other estimates and a linear ephemeris by more than three hours. Note that the vertical scales 
differ between panels. 
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3. Dynamical Models of the Kepler- 11 Planetary System 



Transits of a planet on a keplerian orbit about its star must be strictly periodic. In contrast, 
the gravitational interactions among planets in a multiple planet system cause orbits to speed 



up and slow down 



py small amounts, lead ing to deviations from exact periodicity of transits 



(IHolman & Murray 



2005 



Agoletal 



20051) . Such variations are strongest when planetary orbital 



periods are commen surate or nearly so, which is the case for the large planets Kepler-9 b and c 



(IHolman et al. 



20101) . or when plan ets orbit close to one another, which is the case for the inner 



five transiting planets of Kepler- 1 1 (|Lissauer et al 



2011a). 



To integrate planetary motions, we adopt the 8th order Runge-Kutta Prince-Dormand method, 
which has 9th order errors. Our choice of dynamical epoch was T = 680 days, near the midpoint 
of the fourteen quarters of Kepler data being modeled. In all of our simulations, the orbital period 
and phase of each of the planets are free parameters. The phase is specified by the midpoint of 
the first transit subsequent to our chosen epoch. Initially, we keep all planetary masses as free 
parameters. In some cases, we required planets to be on circular orbits at epoch, whereas in others 
we allowed the orbits to be eccentric. 

We have assumed co-planarity, i.e., negligible mutual inclinations between planetary orbits, 
in all of our dynamical models. We make no attempt to model transit durations or impact 
parameters in our dynamical simulations. 

Our integrations produce an ephemeris of simulated transit times, C s , and we compare these 
simulated times to the observed TTs. We employ the Levenberg-Marquardt algorithm to search 
for a local minimum in x 2 ■ The algorithm evaluates the local slope and curvature of the x 2 surface. 
Once it obtains a minimum, the curvature of the surface is used to evaluate error bars. Other 
parameters are allowed to float when determining the error limits on an individual parameter's 
error bars. Assuming that the x 2 surface is parabolic in the vicinity of its local minimum, its 
contours are concentric ellipses centered at the best-fit value. The orientations of these ellipses 
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depend on correlations between parameters. The errors that we quote account for the increase in 
uncertainty in some dimensions due to such correlations. 

We adopted a wide variety of initial conditions for comparison, and found that our solutions 
were insensitive to the mass of the outer planet, Kepler- 1 1 g. Hence for all subsequent simulations 
used to determine the masses and orbital parameters of the five inner planets, we keep the mass of 
Kepler- 1 1 g as a fixed parameter set to 2.53 x 10 -5 M* (comparable to the masses of similar size 



planet s in this system and equal to 8 M ffi for the value of stellar mass estimated in Lissauer et al. 



201 la|) . with its orbital eccentricity fixed at zero, since eccentrity and mass are inversely correlated 
and our goal is to determine an upper bound on Kepler- 1 1 g's mass. We find that the masses and 
orbital parameters of planets Kepler- 1 1 b-f converge to the values listed in Tables [6] -[8] (Appendix 
B), and the resulting modeled TTVs fit the data well, as displayed in Figures |2]-|7] 

Our dynamical fitting of the planetary parameters minimizes residuals by adjusting 
parameters to search for a best-fit, which is determined by a local minimum value of x 2 - 
Uncertainties are based on the assumption that the shape of the x 2 surface is well- approximated 
by local gradients near the minimum, i.e., is shaped like a parabola. For multi-variate problems 
such as this, the dimensionality of phase space is large, and multiple minima typically exist. 
Furthermore, the low S/N of some lightcurves, particularly, Kepler- 11 b, makes the x 2 surface 
fairly rough, with many local minima. Thus, the minimum that the code finds need not be the 
global minimum, i.e., the best-fit to the data. And even if it does converge to the global minimum, 
parameters that yield other minima with x 2 only slightly larger than that of the global minimum 
are almost as likely to well- approximate the true parameters of the system than are those of the 
global minimum. To qualitatively account for the increased uncertainty caused by these concerns, 
we combined the solutions with the three data sets by averaging their nominal values and defining 
error bars such that they extend over the entire range given by the union of the la confidence 
intervals of all three solutions; error bars are thus asymmetric. Note that this gives fairly large 
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Table 1: Our combined fit dynamical model to the observed transit times, with the orbital periods 
(second column), time of first transit after epoch (third column), ecosu; (fourth column), esinw 
(fifth column), and planetary mass in units of the stellar mass (sixth column), all as free variables 
for planets Kepler- 11 b-f. Periods are given as viewed from the barycenter of our Solar System. 
Because Kepler- 1 1 is moving towards the Solar System at 57 km/s, actual orbital peri ods in the rest 



fram e of Kepler- 11 are a factor of 1.00019 times as long as the values quoted by (|Lissauer et al. 



201 la|) . The simulations used to derive these parameters adopted a circular orbit and a fixed mass 
of 25.3 x 10~ 6 M+ for Kepler- 1 1 g. The upper limit on the mass of planet g was explored separately, 
as described in the text. 
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Fig. 2. — Observed and simulated transit timing variations for planets Kepler- 1 1 b, c and d, using 
transit measurements from E.A. The panels on the left-hand side compare — Q (open symbols 
with error bars) and model times, Q — Q (filled black points). The right hand side plots the resid- 
uals with the dynamical model subtracted from the observed transit timing variations. Note the 
differences between the vertical scales of the various panels. 
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Fig. 3. — Observed and simulated transit timing variations for Kepler- 11 e, f and g, using transit 
time measurements from E.A. See the caption to Figure 2 for details. 
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Fig. 4. — Observed and simulated transit timing variations for Kepler- 11 b, c and d, using transit 
time measurements from J.R. See the caption to Figure 2 for details. 
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Fig. 5. — Observed and simulated transit timing variations for Kepler- 11 e, f and g, using transit 
time measurements from J.R. See the caption to Figure 2 for details. 
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Fig. 6. — Observed and simulated transit timing variations for Kepler- 11 b, c and d, using transit 
time measurements from D.S. See the caption to Figure 2 for details. 
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Fig. 7. — Observed and simulated transit timing variations for Kepler- 11 e, f and g, using transit 
time measurements from D.S. See the caption to Figure 2 for details. 
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ranges, and thus more conservative values than standard la ranges - this is to compensate for 
shortcomings of Levenberg-Marquardt fitting of such a complex multi-parameter space. 

The principal results of our dynamical analysis are presented in Table [U These dynamical 
measurements are combined with estimates of the star's mass and radius to yield measurements 
of the planetary characteristics that we present in Section 5. We also performed fits to each of the 
three sets of TTs with both the eccentricity and the mass of Kepler- 1 1 g were allowed to float, and 
ones in which the mass of planet g was a free parameter but it was constrained to be on a circular 
orbit. In all six cases, the fits converged to values similar to those in our fits with planet g on a 
circular orbit at the nominal mass, albeit with large uncertainties in g's mass, and when allowed to 
float the eccentricity of planet g, and all six of these fits were inferior (in a x 2 /d.o.f. sense) to fits 
with g's parameters fixed. 

To constrain the mass of Kepler- 1 1 g, we performed a suite of simulations using the same 
initial conditions as our best fit to each set of transits times (see Tables [61 [71 [8]). In these fits 
eccentricities for all planets except g were allowed to float, but g's eccentricity is always fixed at 
zero. For each simulation, the mass of planet g was fixed, adding one degree of freedom, to our 
best fit models. The F-ratio, defined as 

A x 2 /A(d.o.f.) 
F_ratl ° = 2 it a n » 

x 2 /(d.o.f.) 

describes the likelihood that a change in the minimum of x 2 could happen by chance given a 
change in the degrees of freedom (d.o.f.), in our case, by fixing the mass of Kepler- 1 1 g to higher 
values. Figure [8] shows the change in \ 2 with variations in the mass of planet g. The 2a limits 
constrain the mass of g, with a confidence of 95%, such that M p (g) < 70 x 10" 6 for two of the 
three datsets. 

One may ask whether as compact a planetary system as Kepler- 1 1 is dynamically stable on 
gigayear timescales. We performed a numerical simulation of a system consisting of planets with 
masses and components of eccentricity la above the nominal values in our best fit (Table [I])- The 
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Fig. 8. — The goodness of fit of our dynamical model to the observed TTs is shown as a function 
of the mass of planet Kepler- 1 1 g. For each point, the x 2 minimum was found keeping the time of 
the first transit after epoch, orbital periods, eccentricities and masses as free variables for planets 
Kepler- 11 b-f. For Kepler- 11 g, the time of its first transit after epoch and its orbital period were 
free parameters, with its eccentricity fixed at zero, and its mass fixed in each numerical run. The 
vertical axis marks the F-ratio, described by Equation (d). Results are shown for the A.E. data 
in green, for the J.R. dataset in blue, and for the D.S. dataset in red. The horizonatal lines mark 
the confidence intervals that \ 2 is not elevated by chance. For the 2a limit, two of the datasets 
constrain the mass of planet g below ~ 70 x lO^M* with 95% confidence. 
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system remained bounded with no gross changes in orbital elements for the entire 3 Myr years 
simulated. 

We performed precise short-term integrations of the nominal system given in Table [T] for 
10 7 days using a Bulirsch-Stoer code. The eccentricities of each of the three low-mass planets, 
Kepler-1 1 b, c and f, varied from minima of ~ 0.002-0.008 to maxima between 0.04 and 0.05. 
The eccentricities of Kepler-1 1 d and e varied from values below 0.0006 to ~ 0.013. Kepler-1 1 g 
was included in these integrations, but it is weakly coupled to the other planets, and its eccentricity 
remained below 0.0006. We also ran an analogous integration with all planetary eccentricities 
initially set to zero. All eccentricities remained small, with peak values for the inner five planets 
in the range 0.0014 - 0.0024. 



4. Properties of the Star Kepler- 11 



1996; 



Lissauer et al. (201 la) performed a standard SME spectroscopic analysis (|Valenti & Piskunov 



Valenti & Fischerll2005|) of a high resolution (R = 60,000) spectrum of Kepler- 11 with 



a wavelength coverage of 360 - 800 nm that was taken by the Keck I telescope a t BJD = 



24555 21.766573 using the observing setup of the California Planet Search group (|Marcy et al. 



20081) . They derived an effective temperature, T eS = 5680 ± 100 K, surface gravity, log g = 4.3 ±0.2 
(cgs), metallicity, [Fe/H] = 0.0 ±0.1 dex, and projected stellar equator ial rotation v sin i = 0.4 ± 0.5 



km/s. Combining these measurements with stellar evolutionary tracks (jGirardi et al 



2000 



Yi et al. 



200 II) yielded estimates of the star's mass, = 0.95 ± 0. 10 M , and radius, = 1 . 1 ± 0. 1 R 



©■ 



We have performed new SME analyses of the same Keck spectrum and of another spectrum 
of comparable quality taken with the same system at BJD = 2455455.802793. The combined 
results (weighted mean values) are: r eff = 5666 ± 60 K, surface gravity, log g = 4.279 ± 0.071 
(cgs), metallicity, [Fe/H] = 0.002 ± 0.040 dex, and projected stellar equatorial rotation 
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vsinz = 3.86 ±0.85 km/s. These values, together with Yale-Yonsei stellar evolutionary tracks, 
yield estimates of the star's mass, = 0.975 ±0.031 M , radius, R+ = 1.193 ±0.1 15 and age = 
9.7 ± 1.5 Gyr. 

The TTV dynamical solution presented in Table [Dprovides stringent constraints on the orbits 
of the inner five transiting planets. We used the computed values of the planets' ecosu and es'mtu 
shown in Table \T\ as constraints in our transit model to provide a geometrical determination of 
the stellar density, p*. The transit model is similar to that described in Appendix A, but we also 
fit for ecosu; and esinui for each of the five inner planets. Posterior distributions for each model 
parameter where estimat ed u s ing a M onte-Carlo-Markov-Chain (MCMC) algorithm similar to the 



one that is described in lFordl (|2005l) . but augmented with a parameter buffer to allow jumps that 
account for correlated variables as described in Rowe et al. (2013, in preparation). We produced 
4 Markov-Chains, each with a length of 400,000. We ignored the first 50% of each chain as burn 
in and combined the remainder into one chain. We adopted the median value for each model 
parameter, which we list in Table |2] 

Since the dynamical model provides a good solution for the orbits of the planets from 
modeling of the TTVs, we reran the transit model and use the constrains on ecosu; and es'mu to 
estimate p*. This translates into a tight constraint on p* = 1.158 ±0.034 g/cm 3 . We combined 
this estimate of p* with the new (weighted mean) SME spectroscopic parameters to determine the 
stellar mass and radius by fitting T e ff, logg and [Fe/H] to M+, age and heavy element mass fraction, 
Z, as provided by the Yale-Yonsei evolution models. We used our MCMC algorithm to determine 
posterior distributions of the stellar parameters and adopted the median value for each parameter 
as listed in Tabled Note that the star is slightly evolved, somewhat more than halfway through its 
lifetime on the main sequence. 

We also conducted a search for spectral evidence of a companion star. We began by fitting 
the observed spectrum of Kepler- 11 obtained on BJD = 2455521.766573 (UT = 21 Nov 2010) 
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Table 2: Transit constraints on the planets of Kepler-11, following dynamical models; b signifies 
impact parameter, i inclination of the orbit to to the plane of the sky and a the orbital semimajor 
axis. 
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Table 3: The characteristics of the star Kepler- 1 1, with la uncertainties. 



with the closest-matching (in a x 2 sense) member of our library of 800 stellar spectra. The stars 
in our library have r eft = 3500-7500 K and logg = 2.0-5.0, which spans the FGK and early M 
type main sequence and subgiant stars. All library stars have accurate parallax measurements, 
allowing for good estimates of stellar mass and radius for each. The Kepler- 1 1 spectrum is placed 
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on a common wavelength scale and normalized in intensity. The x 2 value is then calculated as the 
sum of the squares of the differences between the Kepler- 1 1 spectrum and each library spectrum. 
The final stellar properties are determined by the weighted mean of the ten library spectra with 
the lowest x 2 values. We adopt errors in each parameter by comparing results for standard stars. 
The closest-matching spectrum is modified superficially by removing the Doppler shift relative 
to the observed spectrum, applying needed artificial rotational broadening, setting the continuum 
normalization, and diluting the line strengths (due to a possible secondary star), thereby achieving 
a best-fitting spectrum that can be subtracted from the observed spectrum to yield residuals. 

We search for secondary stars by taking the residuals to that first spectral fit and performing 
the same x 2 search for a "second" spectrum that best fits those residuals; details will be presented 
in Kolbl et al. (in prep). Our approach assumes that spectra are single, until proven double, rather 
than immediately doing a self-consistent two-spectrum fit. This stems from an Occam's razor 
perspective; the notion is that if one spectrum adequately fits the spectrum, without need to invoke 
a second spectrum, then the spectrum can only be deemed single. A minimum in x 2 as a function 
of Doppler shift for the fit of any library spectrum (actually a representative subset of them) to the 
residuals serves to indicate the presence of a second spectrum. We adopt a detection threshold that 
is approximately a 3 a detection of the secondary star. 

We find no stellar companion to Kepler- 1 1 within 0.4" of the primary star, corresponding to 
half of the slit- width (0.87") of the Keck-HIRES spectrometer. The detection threshold for any 
companion star depends on the RV separation between the primary star and the putative secondary 
star. For all RV separations greater than 20 km/s, we would detect (at 3a) any companions that 
are 2% as bright (in the optical) as the primary star. For RV separations of 10 km/s, the detection 
threshold rises to 3% as bright as the primary star, and for RV separations smaller than 10 km/s, 
the detection threshold rises rapidly to unity for FGK stars, but remains at 3% for M dwarfs due 
to their very different spectra. The poor detectability of FGK-type companion stars having little 
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Doppler offset is caused by overlap of the absorption lines. 

Speckle images for Kepler- 11 show no nearby star. Neighbors located in an annulus from 
0.05 to 0.7" from Kepler- 1 1 would have been detected if their brightness were within 3 magnitudes 
in either V or I band, and those between 0.7 and 1.9" distant would have been seen down to a 
magnitude difference of 4 in either band. 



5. Properties of the Planets Orbiting Kepler-11 
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Table 4: The planets of Kepler- 1 1 . The mass and eccentricity of Kepler- 1 1 g are 2a upper bounds. 
All other uncertainties are la confidence intervals. 



Combining the dynamical results presented in Table [Dplus upper bounds on the mass of 
Kepler-1 1 g illustrated in Figure [8]and on this planet's eccentricity from transit models with transit 
parameters of all planets given in Table |2]and the stellar characteristics listed in Table 13 we derive 

d, e a nd 



the planetary parameters shown in Table 4. The nominal mass values of pl anets Kepler- 
f derived herein are within 1 a error bars of the preferred fit presented by Lissauer et al. 



r boiiah . 



and the newl y-estimated masses Keple r-11 b and c are within 2a of their values; the various fits 



presented by 



Migaszewski et al 



the results presented by 



Lissauer et al 



(120 12[) are of comparable accuracy. The major differences from 



(|2011a|) are that the planetary radii are ~ 10% smaller than 
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previously estimated, and planets Kepler- 11 b and especially c are less massive than estimates 
computed with Q1-Q6 data, resulting in the nominal masses monotonically increasing with 
planetary radii rather than the inner pair of planets being more dense than the outer ones. Despite 
the reductions in size estimates, all planets are large for their masses in the sense that they lie 
above both the M p /M m ps (i? p /i? e ) 2 06 relationship that is valid for planets in our Solar System 



(ILissauer et al.ll201 lb|) and mass-radius fits to exoplanets (|Wu & Lithwick 



2013). 



20121 : 



Weiss et al. 



The six planets in Kepler- 11 are all substantially le ss dense than an iron -free ro cky planet, a 



Lissauer et al 



(2011a) and 



Lopez et al 



characteristic already noted for the five inner planets by 
(|2012l) . and which now can be stated with even greater (statistical) significance. As a result, 
they must have substantial envelopes of light components, most likely dominated by the 
cosmically-abundant constituents H 2 , He, and/or H 2 0. In order to understand these envelopes 



we use the thermal evolution models described in detail in Lopez et al. 



( 2012 ). This allows us to 



determine the size of the H/He envelope for each planet, assuming an Earth-like rock/iron core. 



Figure |9] plo ts an updated version of the mass radius diagrams shown in 



(2011a) and 



Lopez et al. 



Lissauer et al 



(120121) . We include all transiting planets with measured masses 
M p < 15 M ffi . For comparison, we include mass-radius curves for Earth-like, 50% water, and 
100% water compositions. In addition, for each of the five Kepler-11 planets whose mass has 
been measured, we include a H/He mass radius curve at the composition and incident flux of that 
planet. 

The new masses imply that Kepler- 1 1 c is less massive than if it were composed of pure 
water, meaning that it must have a large H/He envelope. However, Kepler-11 b can still be 
explained either a H/He or a steam envelope on top of a rocky core. If we assume that Kepler- 1 1 
b's envelope is water rather than H/He, then this planet would be 54%±lt% water by mass. The 
envelope would be composed of steam, since planets like Kepler- 1 1 b are far too irradiated for 
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their interiors to include liquid or high-pressure ice phases. Most of the H 2 would be in the vapor 
and molecular fluid phase s, with the ionic fluid and pla sma phases occurring at high pressures 



deep within these planets (INettelmann et al.1 



2008 



2011b 



For mixtures of rock with H/He (no H 2 0), and using the sizes and masses presented in 
Table 4, we find that Kepler- 1 1 b is currently 0.41%±g;g| H/He, Kepler- 1 1 c is 4.6%±^ H/He, 
Kepler- 11 d is 6.0%±} H/He, Kepler- 11 e is 13.9%±}J| H/He, and Kepler- 11 f is 3.5%±?;^ 
H/He by mass. The quoted uncertainties include the measured uncertainties on each planet's 
mass, radius, incident flux, and age as well as theoretical uncertainties on the albedo and the iron 



fraction of the rocky /iron core ([Marcus et al 



2010)) . Despite the small mass fractions in light 



gases, the presence of these H/He envelopes is key to the observed radii. One way to emphasize 
this fact is to compare each planet's radius to the radius of its rocky core, as determined by our 
thermal evolution models for planets lacking H 2 0. For every Kepler- 1 1 planet whose mass has 
been measured except for b, approximately half of the observed radius is due to its H/He envelope. 
The cores make up 47%, 54%, 42%, and 50% of the total radii of planets Kepler-1 1 c, d, e, & f, 
respectively, and thus only 7 - 16% of the volume. Moreover, even for Kepler-1 1 b, the rocky core 
only makes up 67% of the total radius, corresponding to 30% of this planet's volume. 

In addition, we h ave included a updated version of the mass loss threshold diagram presented 



in 



(|2012l) . Figure [TOl plots incident flux against the product of planet mass times 
planet density. Lines in this space correspond to mass loss timescales, making this diagram 
useful for underst anding how the population of highly-irradiated planets has been sculpted by 



photoevaporation (|Lecavelier des Etangs 



|2007|) . Here we have color-coded planets by the fraction 
of their mass in the H/He envelope, assuming an Earth-like core. Four known exoplanets are dense 
enough to be composed of bare rock; these planets are shown as rust colored. The key feature 
of Figure [TOl is that there is a critical mass loss timescale above which there are no planets with 
significant H/He envelopes. The dashed black line shows the critical mass loss timescale found by 
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(ILopez et al.ll2012|) . The existence of such a mas s loss threshold is a robust predication of p lanet 



evolution models that include photoevaporation (|Owen & Jackson 



20121 : 



Lopez et al 



2012). The 



three planet s that lie above th i s threshold in the upper right are 



CoRo T-7 b (ILeser et al 



2011 



Demory et al 



2009 



Oueloz et al. 



2009 



Hatzes et al. 



Kep 



er-lOb (Batalha et al. 



201 1), and 55 Cancri e (Win n et al 



2011) 



201 lh . none of which are expected to have H/He envelopes. 



Lopez et al. 



(120121) 



With the new masses, Kepler- 11 b & c now lie slightly above the 
threshold. This does not affect the general conclusion that photoevaporation has sculpted the 
population of super-Earths and sub-Neptunes found at short period. A single critical mass loss 
timescale is only a rough approximation. In reality mass loss evolution depends on many unknown 
parameters like the size of initial accreted H/He envelope, which can easily produce a factor of ~ 
2-3 scatter in the mass loss threshold (Lopez & Fortney 2013, submitted). However, it curious 
that planets Kepler- 1 1 d, e, & f all lie on the same mass loss timescale, while planets Kepler- 1 1 
b & c seem to lie on a different timescale. This suggests that either Kepler- 1 1 b & c formed with 
larger initial envelopes or that they experienced mass loss with significantly lower efficiency. 
The latter case is a definite possibility, since at early times planets Kepler- 11 b & c experienced 
EUV fluxes that put them rig ht in the transition betwee n energy-limited and the less efficient 



recombination-limited escape (|Murray-Clay et al. 



2009) 



6. Conclusions 

We have performed an updated analysis of the Kepler- 1 1 planetary system, concentrating on 
the dynamical interactions evident in transit timing variations observed in the first three years of 
Kepler photometric data. We have also improved our estimates of the characteristics of the star 
using new techniques to analyze a spectrum taken at the Keck observatory. Our updated transit, 
stellar and planetary parameters are presented in Tables [3] and 4, respectively. 
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The six planets observed to transit Kepler- 11 all have small orbital eccentricities. None is 
dense enough to be composed entirely of rocky material, and at least the four middle planets must 
contain volumetrically- significant envelopes of gases less dense than H 2 0. Planets Kepler-11 b 
and f, and nominally c as well, are less massive than any other exoplanets for which both mass and 
radius have been measured. The planetary parameters are consistent with a monotonic increase in 
mass as a function of radius, although the Kepler- 1 1 planets are less dense for a given radius than 
most other planets with mass and radius measurements. 

Kepler was competitively selected as the tenth Discovery mission. Funding for this mission 
is provided by NASA's Science Mission Directorate. E. A.'s work was supported by NSF Career 
grant AST-0645416. W. F. W. gratefully acknowledges support from the Kepler Participating 
Scientist Program via NASA grant NNX12AD23G, and from the NSF via grant AST- 1 109928. 
We thank Jerome Orosz and Gur Windmiller for assistance in developing D.S.'s method for 
measuring transit times and Darin Ragozzine and Billy Quarles for helpful comments on the 
manuscript. 

7. Appendix A: Techniques used to Measure Transit Times 

We measured transit times using three different techniques, each of which is described below. 

7.1. TT Measurements by Jason Rowe 

This analysis used Ql - Q14 long-cadence and Q3 - Q14 short cadence Kepler simple 
aperture photometry (labeled SAP_FLUX). Only data with a quality flag set to zero as documented 
in the Kepler data release notes were used. This provided 52,539 and 1,464,980 long and short 
cadence photometric measurements, respectively. 
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The data were initially detrended using a running 2-day box-car median filter that was applied 
to individual segments of time-series photometry. A segment was defined as a continuous string 
of time-series data that does not contain an interruption longer than 2.5 hours (5 long cadence 
measurements). This was done to handle offsets observed after data outages, typically due to a 

change in the thermal environment of the CCD detector. A circular quadratic transit model based 

i — ° 1 i 1 1 



on 



Mandel & Agoll (|2002l) was fit to the data by minimization of x 2 with a Levenberg-Marquardt 



algorithm. The transit model was used to measure the transit duration for each transiting 
planet. The original SAP_FLUX photometric data were then reprocessed using a second-order 
polynomial to detrend the time-series to remove instrumental (such as focus changes) and 
astrophysical effects. All data obtained during transit were excluded, as well as those taken in the 
30 minutes before ingress and in the 30 minutes after egress. A clipping algorithm was used to 
exclude any measurement that differed from the mean by more than 3a. Measurements obtained 
during a planet transit were excluded from the clipping exercise. It was found that the data before 
a data outage near JD = 2455593 could not be sufficiently detrended. As such, data from 2455593 
to 2455594.5 were excluded, which meant that a transit of Kepler- 1 1 g was not included in our 
analysis. 

The detrended LC and SC photometric time-series were then each fit with a multi-planet, 
circular orbit, quadratic Mandel & Agol transit model. The model parameters are the mean stellar 
density (p*), photometric zero point, and, for each planet, the center of transit time, orbital period, 
impact parameter, and scaled planetary radius (R p /R+). The model assumes that the mass of star is 
much greater than the combined mass of the orbiting planets, so that 



a \ it (M^+M p ) 



RJ 3GP 2 

5 * 



P*. (2) 



A photometric time- series for each transiting planet was then produced by removing the 
transits of the other transiting planets. The remaining transits were then individually fit by using 
the best-fit model as a template and only allowing the center of transit time to vary. This yielded 
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a time-series of transit timing variations (TTVs) for each planet. The measured TTVs where 
then used to linearize (or deTTV) the photometry, such that when folded at the orbital period the 
transits are aligned in the resulting lightcurve. The multi-planet transit model was then refit out to 
the deTTVed lightcurve and used the updated template to determine the final set of TTVs shown 
by the green points in Figure 1 . Uncertainties in the transit times were determined by examining 
the residuals from the fits to each individual transits and scaling the photometric errors such that 
reduced x 2 wa s equal to one. The diagonal elements of the co-variance matrix were adopted as 
the uncertainty in the measurement. 

7.2. TT Measurements by Eric Agol 

The times of transit were fit using a quadratic limb-darkening model in which the duration 
and impact parameter for each planet were assumed to be fixed, while the times of each transit 
were allowed to vary. The model was computed simultaneously for the short cadence (when 
available) and long cadence data (otherwise). A window equal to one transit duration was included 
before and after every transit. The lightcurve was divided by the model (computed for all planets 
simultaneously so that overlapping transits were properly accounted for), and then fit with a 
third-order polynomial for each contiguous data set (without gaps larger than 12 hours). The 
model parameters were optimized until a best fit was found; a second iteration was carried out 
after outliers from the first fit were rejected. After finding the best fit, the times of each and every 
transit were allowed to vary over a grid of values spanning (typically) about 2 hours on either side 
of the best fit time. The variation in x 2 with transit time was then fit with a quadratic function 
to measure the uncertainty in the transit time. If that fit failed, then the transit time error was 
measured from the width of the x 2 function for values less than one above the best fit value. 
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7.3. TT Measurements by Donald Short 

In contrast to the Rowe and Agol methods, a purely mathematical technique was used to 
determine the transit times, under the assertion that the time of a transit event can be estimated 
without need of a physical model of the event. Under conditions of poor signal-to-noise ratio 
or undersampling, the constraints imposed by a physical model are extremely valuable. For 
high signal-to-noise cases, a non-physical model can match, or even excel a physical model 
under certain conditions. The limitations in a physical model, such as imperfect limb darkening 
parameterization or assumed zero eccentricity, have no consequence in a non-physical model. 
Since no assumptions about sphericity, obliquity, gravity darkening, strict keplerian motion, etc., 
were made, the method is insensitive to errors in these physical parameters or effects. 

Both LC and SC data were used in computing the planetary transit time estimates, provided 
the pipeline data quality flag had the nominal value of zero. The TTs were estimated by an 
iterative method starting with the SC data. Using an estimate of the transit du ration and estimates 



Lissauer et al. 



(|2011al) . each transit was 



of the transit times based on the linear ephemeris from 
locally detrended. Detrending employed a low-order polynomial centered on the transit and 
extending symmetrically either 0.3, 0.6, or 0.83 days beyond the ends of the transit; the length and 
polynomial order that provided the best fit to these out-of-transit data was selected. During this 
process, each transit was checked for missing data and overlapping transits from other planets that 
could compromise the determination of that TT. Transits that had such problems were eliminated 
from further consideration. After detrending, the transits were shifted in time so that the center 
of each transit was at time zero. All of the transits were then combined ("stacked" or "folded" 
on top of each other). A piecewise cubic Hermite spline (PCHS) was then least-squares fit to 
the combined-transit lightcurve, giving a transit template. The transit template was generated by 
the data themselves; no physical constraints on its shape were imposed. As such, it should be an 
excellent match to the observed transits. From this template, a refined transit width was estimated 
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and used to revise the detrending of each transit. The template was then correlated with each 
individual transit, yielding improved TTs. Any outliers with respect to the template were flagged 
and eliminated from further template building, but no rejections were made when estimating the 
individual TTs. The detrended transits were shifted (folded) on the revised TTs, combined, and 
a new PCHS template generated. Again, the individual transits were then detrended, now using 
both the revised duration and revised transit times. The detrended transits were correlated with the 
revised template, yielding a refined set of TTs. Three iterations of this process were carried out. 
The uncertainty in each TT was estimated from the shifts in time needed to degrade the x 2 fit of 
the template to the transit by one. 

For transits with LC data only, the SC PCHS template was convolved to 30 minutes, yielding 
the LC template. The LC template was then correlated with each transit, providing a correction to 
the times from the initial linear ephemeris. The revised TTs were used to improve the detrending 
window, but the template was not updated - it was held fixed at the shape derived from the SC 
template. This process iteratively produced measurements of the TTs, uncertainties, and model 
fits for each transit. Finally, those TTs that had large timing error estimates (>40 min) were 
eliminated from the final list of TTs. 

The process above was repeated independently for each planet, noting that overlapping 
transits from different planets were discarded. In general, the TTs computed by this method agree 
quite well with the physical methods; however, the error estimates are notably larger. 

8. Appendix B: Details of Dynamical Models 

Here we present the results of our dynamical models in detail. We carried out three classes 
of fit using each set of TTs. In the "all-circular" class, all planets were assumed to travel on 
circular orbits at epoch. In the "all-eccentric" class, all planets were allowed to have eccentric 
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orbits at epoch. We found that quality of these fits were not sensitive to the mass or eccentricity 
of planet Kepler- 1 1 g as long as these were not too large, so we performed "g-fixed" fits wherein 
the eccentricity of planet g is set to zero at epoch and its mas s set to 25.3 x 10~ 6 M + (which equals 



8 A/ e for assume a stellar mass of 0.95 Af©, as estimated by 



Lissauer et al. 



2011 



Table [5J compares the quality of fit between using various data sets and assumptions. Note 
that comparisons of the numerical values between the quality of fits using different data sets are 
not meaningful because of the differing prescriptions employed to compute the uncertainties of 
individual TTs, but comparison between the reduced x 2 f° r me all-circular, all-eccentric, and 
g-fixed results using a given set of TTs shows that eccentricities are detected for the five inner 
planets but not for planet g. As the quality of the all-circular fits are distinctly inferior to those that 
allow at least the five inner planets to travel on eccentric orbits, we do not consider the all-circular 
fits further. 

As shown in Table [5j the g-fixed fits, which are presented in Tables [61 [7] and [81 are of 
slightly better quality (in a x 2 /d.o.f. sense) than are the corresponding all-eccentric fits. Thus, the 
parameters from the three g-fixed fits are synthesized to incorporate the full ranges of all 1 a error 
bars from fits to each set of data and displayed as our primary results in Table Q] Table [9] is the 
counterpart of Table 1, synthesizing all-eccentric fit results of the three sets of transit time data. 

The small values (compared to unity) of x 2 /(d.o.f.) shown in Table[5jfor fits to D.S.'s TTs 
imply that the uncertainties quoted for these TTs were overestimated. Similarly, the large values 
of x 2 /(d.o.f.) for E.A.'s TTs strongly suggest that these uncertainties were underestimated. The 
values of % 2 /(d.o.f.) near unity for both fits allowing eccentric planetary orbits to J.R.'s TTs 
suggest that uncertainties in these TTs may have been slightly overestimated. 
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Planet 


All 


circular 




All 


eccentric 




g- fixed 




(Tables HEM 




Xea 


Xm 


Xds 


Xea 


Xm 


Xds 


Xea 


X 2 m 


Xds 


b 


229.62 


100.36 


41.61 


189.49 


86.59 


36.16 


189.43 


86.61 


36.14 


c 


280.39 


111.72 


72.30 


211.18 


79.71 


50.15 


211.21 


79.66 


50.14 


d 


123.25 


51.68 


15.63 


100.95 


45.08 


14.78 


101.02 


45.08 


14.79 


e 


78.47 


53.97 


22.06 


46.47 


29.08 


12.03 


46.73 


29.24 


12.06 


f 


120.79 


58.09 


38.70 


52.03 


16.81 


12.73 


54.10 


17.24 


13.06 


g 


9.84 


6.71 


3.69 


10.08 


6.69 


3.62 


9.61 


6.66 


3.64 


total 


842.16 


382.53 


193.99 


610.20 


263.97 


129.46 


612.09 


264.49 


129.83 


X 2 /(d.o.f.) 


3.25 


1.38 


0.82 


2.47 


1.00 


0.58 


2.46 


0.99 


0.57 



Table 5: x 2 contributions from each planet for a suite of models against both sets of transit times. 
The second through fourth columns show best fits to an orbital configuration with all eccentricities 
fixed at zero, the fifth through seventh columns show all eccentric fits, and the eighth through tenth 
columns shows the g-fixed models whose results are shown in Tables |6l [7] 

and [8] 



Planet 


P (days) 


To 


ecosa> 


es'mco 


M p /M± x 10" 6 


b 


10.3043±0.0002 


689.7378±0.0009 


0.038 ±0.016 


0.009 ±0.008 


3.91±1.03 


c 


13.0236±0.0003 


683.3494±0.0005 


0.019±0.014 


-0.005 ±0.004 


6.23±1.75 


d 


22.6839±0.0003 


694.0061±0.0005 


-0.006 ±0.003 


0.001 ±0.001 


23.60±1.66 


e 


31.9996±0.0004 


695.0752±0.0005 


-0.009 ±0.002 


-0.009 ±0.001 


27.77±1.92 


f 


46.6903±0.0011 


718.2737±0.0015 


0.007 ±0.003 


-0.007 ±0.002 


7.45±1.09 


g 


118.3807±0.0004 


693.8022±0.0010 


(0) 


(0) 


(25.29) 



Table 6: Best (g-fixed) dynamical fit to TTs from E.A. 
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Planet 


P (days) 


T 


e cos 00 


esinu 


M p /M+ x 10" 6 


b 


10.3039±0.0004 


689.7391±0.0012 


0.050 ±0.019 


0.014±0.010 


6.80±2.16 


c 


13.0240±0.0005 


683.3497±0.0010 


0.033 ±0.016 


0.005 ±0.008 


9.25±3.34 


d 


22.6849±0.0005 


694.0072±0.0007 


-0.003 ±0.004 


0.004 ±0.003 


23.61±1.84 


e 


31.9999±0.0005 


695.0756±0.0005 


-0.007 ±0.003 


-0.007 ±0.003 


23.85±2.55 


f 


46.6877±0.0014 


718.2697±0.0021 


0.014 ±0.005 


-0.001 ±0.002 


5.26±1.21 


g 


118.3806±0.0005 


693.8010±0.0010 


(0) 


(0) 


(25.29) 



Table 7: Best (g-fixed) dynamical fit to TTs from J.R. 



Planet 


P (days) 


To 


e cos uj 


esinw 


M p /M± x 10" 6 


b 


10.3036±0.0007 


689.7363±0.0032 


0.009 ±0.008 


0.072 ±0.018 


6.80±3.28 


c 


13.0247±0.0006 


683.3490±0.0015 


-0.004 ±0.005 


0.059 ±0.014 


12.06±6.25 


d 


22.6846±0.0006 


694.0074±0.0016 


-0.001 ±0.002 


-0.000 ±0.001 


21.27±3.23 


e 


31.9993±0.0009 


695.0759±0.0011 


-0.008 ±0.003 


-0.010 ±0.004 


22.91±4.73 


f 


46.6883±0.0027 


718.2695±0.0023 


0.014 ±0.007 


-0.007 ±0.005 


5.94±2.55 




118.3809±0.0003 


693.8030±0.0021 


(0) 


(0) 


(25.29) 



Table 8: Best (g-fixed) dynamical fit to TTs from D.S. 
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Planet P (days) T ecosu esinu M p /M+ x 10 

b 10 3039+ 0006 689 7377+0 0031 n 032+ 037 032+ 060 5 S3+ 4 - 29 

U 1U.OUO^_ 00()11 OO^./O/ /_oo046 u - u:5Z -0.035 u - u:)z -0.030 *' ,o,; '-3.09 

c 13.0241S 683.3494^ 0.016£8i 0.02Q™* 9.133$ 

d 22.6845^;^ 694.0069^ -0.003^ 0.002™« 22.84^;^ 

e 31.9996^ 695.0755^ -O.OMflffi -0.009$% 24.83^ 

f 46.6887^^1 718.2711^ 0.01 1^ -0.005^ 7 6 6.20*;* 

g U8380n8§tg 693.8021t°°g° 0.032$°% 0.022$°% 232Vjjjl 



Table 9: Dynamical all-eccentric fits to the observed transit times with the orbital periods (second 
column), time of first transit after epoch (third column), ecosa; (fourth column; e represents eccen- 
tricity, co is the angle, measured from the star, between the place the planet's orbit pierces the sky, 
coming towards the observer, and the pericenter of the orbit), esinw (fifth column), and planetary 
mass in units of the stellar mass (sixth column), all as free variables. For 'g', this model has settled 
on a mass near the initial estimate of 8M ffi (25.3 x 10" 6 M*). 
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Fig. 9. — Updated mass-radius diagram for transiting exoplanets with measured masses, along 
with curves for different compositions. Planets are color-coded by the incident bolometric flux 
they receive. Kepler planets are shown by circles, filled Kepler- 1 1, open for others, with numbers 
and letters indicating each planet. Other known exoplanets in this mass and radius range are shown 
by open squares; in order of increasing radius, these are CoRoT-7 b, 55 Cancre e, GJ 1214 b and 
GJ 3470 b. Solar System planets Venus and Uranus are shown by black letters. The solid black 



curve is for an Earth-like composition wit 



l 2/3 r ock and 1/3 iron by mass. All other curves use 



thermal evolution calculations (|Lopez et al 



20121), assuming a volatile envelope atop a core of rock 



and iron with composition the same as that of the bulk Earth. The dashed blue curves are for 50% 
and 100% water by mass. The dotted orange curves are for H/He envelopes at 8 Gyr; each one is 
tailored to match a Kepler- 1 1 planet and is computed at the appropriate flux for that planet. 
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Fig. 10. — Updated version of the mass loss threshold diagram from ; 
ric flux at the top of the atmosphere, relative to the flux incident on Earth, is plotted the product 
of planet mass and planet density. Again, the Kepler-11 planets are shown by filled circles with 
letters indicating each planet. Open squares show the other extrasolar planets < 15M®, while 
crosses show all other transiting planets with measured masses up to 100 M©. Planets are color 
coded by the percentage of their mass in their H/He envelopes, according to thermal evolution 
models. Potentially rock y planets are rust colored. The dashed black line shows the critical mass 



loss timescale found by 



Lopez et al. 



d2012h . 
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